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24 CHAPTER 2 Mathematical Models of Systems 

Starting at the output we obtain 

I(s) = G 1 (s)G 2 (s)E(s). 

But E(s) = R(s ) — H(s)I(s), so 

I(s) = Gi(«)G 2 («) [R(s) - H(s)I(s)} . 

Solving for I (s) yields the closed-loop transfer function 

I(s) G 1 (s)G 2 (s) 

R(s) 1 + Gi {s)G 2 (s)H{s) ' 

E2.8 The block diagram is shown in Figure E2.8. 




FIGURE E2.8 

Block diagram model. 

Starting at the output we obtain 

Y(s) = - Z(s ) = -G 2 (s)A( S ). 
s s 

But j4(s) = Gi(s) [— H 2 (s)Z(s) — H^(s)A(s) + W(s)] and Z(s) = sY(s), 
so 

y( s ) = -G 1 (s)G 2 (s)H 2 (s)Y(s) - Gi(s)H 3 (s)Y (s) + -G 1 (a)G 2 (a)W(a). 

s 

Substituting W(s) = KE(s) — Hi(s)Z(s) into the above equation yields 

y(s) = —G\{s)g 2 {s)h 2 {s)y{s) - Gi( s )F 3 ( s )y(s) 

+ Igi(s)G 2 (s) [KE{s) - Hi(s)Z(s)\ 
s 
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and with E(s) = R(s) — Y(s ) and Z(s) = sY(s ) this reduces to 
Y(s) = [-Gi (s)G 2 (s) (. H 2 (s ) + H^s)) - Gi(a)fr 3 (a) 
- ^G 1 (8)G 2 (s)K\Y(s) + ^G 1 (s)G 2 (s)KR(s). 

Solving for Y ( s ) yields the transfer function 

Y(s) = T(s)R(s), 



where 

T(s) = 



KG 1 (s)G 2 (s)/s 



1 + Gi(a)G 2 (a) [(H 2 {s) + H^s)} + G^H^s) + KG^G^s) / s' 
From Figure E2.9, we observe that 

F f (s) = G 2 (s)U(s) 

and 



FrW 

Then, solving for U(s) yields 



G 3 (s)U(s) . 



u ^ = 7W) F ' {s) 

and it follows that 

F * (s) = Wj u(s) ■ 

Again, considering the block diagram in Figure E2.9 we determine 
F f (s) = Gi(s)G 2 (s)[i*(s) - H 2 (s)F f (s) - H 2 (s)F R (s )] . 

But, from the previous result, we substitute for F R (s) resulting in 
F f (s) = G 1 (s)G 2 (s)R(s)-G 1 (s)G 2 (s)H 2 (s)F f (s)-G 1 (s)H 2 (s)G 3 (s)F f (s) . 
Solving for Ff(s) yields 



Ff(s) 



' Gi(a)G 2 (a) ' 

.1 + Gi(s)G 2 (s)FT 2 ( S ) + G 1 (s)G 3 (s)H 2 (s)_ 



R(s) . 
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Problems 51 

P2.28 The signal flow graph is shown in Figure P2.28. 

-d 




FIGURE P2.28 

Signal flow graph. 



(a) The PGBDP loop gain is equal to -abed. This is a negative transmis- 
sion since the population produces garbage which increases bacteria 
and leads to diseases, thus reducing the population. 

(b) The PMCP loop gain is equal to +efg. This is a positive transmis- 
sion since the population leads to modernization which encourages 
immigration, thus increasing the population. 

(c) The PMSDP loop gain is equal to +ehkd. This is a positive trans- 
mission since the population leads to modernization and an increase 
in sanitation facilities which reduces diseases, thus reducing the rate 
of decreasing population. 

(d) The PMSBDP loop gain is equal to -fehmed. This is a positive 
transmission by similar argument as in (3). 

P2.29 Assume the motor torque is proportional to the input current 



Then, the equation of motion of the beam is 

J(f) = ki , 

where J is the moment of inertia of the beam and shaft (neglecting the 
inertia of the ball). We assume that forces acting on the ball are due to 
gravity and friction. Hence, the motion of the ball is described by 

rnx = mgef) — bx 
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Prom the block diagram we have 

n(s) = G 2 (s)[G 1 (s)E 1 ( S ) + G 3 (s)E 2 (s )] 

= G 2 (s)G 1 (s)[R 1 (s) - P^i(s)y 1 (s)] + G 2 (s)G 3 (s)E 2 (s ) . 

Therefore, 

GiOQGaOQ N G 2 (s)G 3 (s) 

l( } 1 + Gi (s)G 2 (s)H 1 (s) Rl{ )+ l + G^G^H^s) ] ' 

And, computing E 2 (s) (with R 2 (s ) = 0) we find 

-G 4 (s 



E 2 {s ) = H 2 (s)Y 2 (s) = H 2 (s)G 6 (s ) 



or 



E 2 (s) = 



[G 2 (s) 
G 4 (s)G 6 (s)H 2 (s) 



Yi(s) + G 5 (s)E 2 (s) 



Yi(s) . 



Yi(s) 



G 2 (s)(1 - G 5 (s)G 6 (s)H 2 (s)) 

Substituting E 2 {s) into equation for Yi(s) yields 

Yl (s) = Gl(g)G2(5) RUs) 

lU 1 + G 1 (s)G 2 (s)H 1 (s) lU 

G 3 (s)G 4 (s)G 6 (s)H 2 (s) 

(1 + G^G^H^il - G 5 (s)G 6 (s)H 2 (s)) 

Finally, solving for Ti(s) yields 

Y 1 {s) = T 1 (s)R 1 (s) 

where 
Ti(s) = 

r Gi( s )G 2 ( s )(i-G 5 ( s )G 6 (s)F 2 (s)) 

.(1 + Gi(s)G 2 (s)Hi(s))(1 - G 5 (s)G 6 (s)H 2 (s)) - G 3 (s)G 4 (s)G 6 (s)H 2 (s) 

Similarly, for T 2 (s) we obtain 

Y 2 (s) = T 2 (s)Ri(s) . 

where 
T*(s) = 

Gi(a)G 4 (8)G 6 (g) 

(1 + Gi(s)G 2 ( S )Pf 1 (s))(l - Gs(s)G 6 (s)H 2 (s)) - G 3 (s)G 4 (s)G 6 (s)H 2 (s) J 
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P2.33 The signal flow graph shows three loops: 



L\ = —G1G3G4H2 
L2 = —G2G5GQH1 
L% = — H1G8GQG2G7G4H2G1 . 

The transfer function Y 2 /R 1 is found to be 

Y 2 (s) GiGgGeAi - G 2 G 5 G 6 A 2 
R\ (s) 1 — (L\ + L2 + T3) + (L\L 2 ) 

where for path 1 

Ai = 1 



and for path 2 

A 2 = 1 — L\ . 

Since we want Y 2 to be independent of R,\ , we need Y 2 /R 1 = 0. Therefore, 
we require 

GrGgGe - G 2 G 5 G 6 ( 1 + G^G^) = 0 . 

P2.34 The closed- loop transfer function is 

Y(s) _ G 3 (s)G 1 (s)(G 2 (s) + K 5 K 6 ) 

R(s) 1 — Gs(s)(Hi(s) + Kq) + G 3 (s)Gi(s)(G 2 ( S ) + K§K§)(H 2 (s) + K±) 

P2.35 The equations of motion are 

mijji + b(y i - y 2 ) + h(yi - y 2 ) = 0 
m 2 m + Hm - yi) + k 1 (y 2 - yi) + k 2 y 2 = k 2 x 

Taking the Laplace transform yields 

(■ m\s 2 + bs + fci)Yi(s) - ( bs + ki)Y 2 (s) = 0 
(m 2 s 2 + bs + k\ + k 2 )Y 2 {s) - (bs + ki)Yi(s) = k 2 X(s) 

Therefore, after solving for Yi(s)/X(s), we have 

l 2 (s) _ k 2 (bs + ki) 

X(s) ( m\s 2 + bs + ki)(rri 2 S 2 + bs + k\ + k 2 ) — (bs + k\) 2 

P2.36 (a) We can redraw the block diagram as shown in Figure P2.36. Then, 

= Ki/s(s + l) = K\ 

1 + ATi(l + K 2 s)/s(s + 1) s 2 + ( 1 + K 2 K 1 )s + K 2 ' 
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Exercises 



E3.1 One possible set of state variables is 

(a) the current il 2 through L 2 , 

(b) the voltage vc 2 across C 2 , and 

(c) the current ij jX through L\. 

We can also choose ucr , the voltage across C\ as the third state variable, 
in place of the current through L\. 

E3.2 We know that the velocity is the derivative of the position, therefore we 
have 

dy _ 
dt V ’ 

and from the problem statement 



dv 

dt 



—k\v(t) - k 2 y(t ) + k 3 i(t) . 



This can be written in matrix form as 




Define u = i, and let k\ = k 2 = 1- Then, 



x = Ax + Bu 



where 



0 1 

-1 -1 



0 


, and x = 


( V \ 


_ ^3 _ 




\ v ) 
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CHAPTER 3 State Variable Models 



E3.3 The charactersitic roots, denoted by A, are the solutions of det(AI — A) = 
0. For this problem we have 



det(AI — A) = det 



A -1 
1 (AH- 1) 

Therefore, the characteristic roots are 



= A(A + !) + ! = A 2 + A + 1 = 0. 



x 1.^/3 , x 1 ->/3 

Ai = --+j — and A 2 = ---j — 



2 " 2 

E3.4 The system in phase variable form is 



x = Ax + B u 

y = Cx 



where 



A = 



0 


1 


0 




0 


0 


0 


1 


, B = 


0 


-8 


-6 


-4 




20 



C = 



1 0 0 



E3.5 From the block diagram we determine that the state equations are 

22 = —(fk + d)x i + ax i + fu 
X\ = — kx2 + u 



and the output equation is 



Therefore, 



V = bx 2 



x = Ax + B'/j 
y = Cx + Du 



A = 



where 

0 —k 

a — {fk + d ) 

E3.6 (a) The state transition matrix is 

At 



, B = 


1 


O 

II 


0 b 




. / . 




L J 



and D = [0] 



&(t) = e At = I + A t+ -A 2 t 2 + 
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In state variable form we have 



x = 



(fc+fci) 

m 



ki 

m 



b_ 

m 



o 

o 



0 

ki 



(fc+fcl) 

m 



where x\ = x, X2 = x, X3 = q and X4 = q. 
E 3.16 The governing equations of motion are 



0 

0 

1 

~b_ 

m 



mix + k\(x — q) + — q) = u(t ) 

m 2 q + k 2 q + b 2 q + b\ (q - x) + k\ (q - x) = 0 . 



Let xi = x, 


X 2 = X 


X 3 = 


q and X4 = 


q. Then, 






0 


1 


0 


0 




0 




ki 


bi 


ki 


bi 




1 


X = 


mi 


m± 


mi 


mi 


X + 


mi 




0 


0 


0 


1 




0 




ki_ 


bi 


(fci+fca) 


(61+62) 




0 




m2 


m2 


m 2 


m2 







Since the output is y(t) = q(t), then 



V = 



0 0 10 



x . 



E 3.17 At node 1 we have 



Civ 1 



Vg - V\ V2 - Vl 

R\ R 2 



and at node 2 we have 



C 2 v 2 



Vb~V 2 vi - v 2 

R3 R2 



Let 



Xl = Vl 



and 



x 2 = v 2 . 
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Then, in matrix form we have 



x = 



( RiCi + RaC’J 

1 



1 



1 



H 2 C i 

+ 



X+ 



R 2 C 2 \R 3 C 2 ' R 2 C 2 

E3.18 The governing equations of motion are 



T di 1 

Rii + L i— + v = v a 
at 

t di 2 , 

L 2-J- + v = V b 
at 

„dv 

lL=l 1 + *2 = C— . 

dt 



1 



Hi Cl 

0 



0 

1 



R 3 C 2 





V a 




. Vb . 



= *2, 


CO 

II 

C5 


Ul 


= v a and U 2 = 


= V b . 


Then, 




R 


0 


1 




1 


0 




Li 


Hi 




Li 


X = 


0 


0 


_ j_ 

L 2 


x + 


0 


j_ 

L 2 




1 

C 


1 

c 


0 




0 


0 


y = 


0 0 


1 


x+[0] 


u . 







u 



E3.19 First, compute the matrix 



si- A = 



s —1 
3 s + 4 



Then, <h(s) is 



4>(s) = (si - A)” 1 = 



A(s) 



s + 4 1 
-3 s 



where A(s) = s 2 + 4s + 3, and 
G(s) = 



10 0 



s+4 1 

A(s) A(s) 

3 s 





1 

O 

1 




1 



10 



s 2 + 4s + 3 



E3.20 The linearized equation can be derived from the observation that sin 8 ~ 6 
when 9 ~ 0. In this case, the linearized equations are 

6 + yd + —0 = 0 . 

L m 
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m 



FIGURE P3.5 

Block diagram model. 




■m 



P3.6 The node equations are 

dv-\ Vj — vi 

0.00025—— +i L - — A = 0 

dt 4000 

dv 9 Vo 

0.0005— - 1 -i L + -i 3 = 0 

dt 1000 

dir 

0 . 002 — — |- V 2 — v\ = 0 
dt 

Define the state variables 

X\ = V\ X2 = v 2 x 3 =il ■ 

Then, 

x = Ax + Bu 



where 



-1 


0 


-4000 




1 


0 


0 


-2 


2000 


, B = 


0 


2000 


500 


-500 


0 




0 


0 



P3.7 Given K = 1, we have 



KG(s) • - 
s 



(s+ l) 2 
s(s 2 + 1) 



We then compute the closed-loop transfer function as 



T(s) 



s 2 + 2s + 1 
3s 3 + 5s 2 + 5s + 1 



s 1 + 2s 2 + s 3 
S + Ss-i + Ss - 2 -^^- 3 ‘ 
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(c) The characteristic equation is 



det [si — A] = det 



s -1 0 

Os -1 

5 1 s + 10.1 



s 3 + 10.1s 2 + s + 5 = 0 . 



The roots of the characteristic equation are 

si = -10.05 and s 2 , 3 = -0.0250 ± 0.7049j . 

All roots lie in the left hand-plane, therefore, the system is stable. 

P3.10 (a) From the signal flow diagram, we determine that a state-space model 

is given by 





' -K] 


k 2 




’ Ki 


K 2 ' 




n 


X = 


~K\ 


-k 2 


x + 


K\ 


K 2 




r 2 



yi 




1 


0 


2/2 




0 


1 



(b) The characteristic equation is 

det [si - A] = s 2 + (K 2 + K\)s + 2AiA 2 = 0 . 

(c) When K\ = K 2 = 1, then 



A = 



-1 1 
-1 -1 



The state transition matrix associated with A is 
$ = £" 1 {[sl- A]” 1 } 



= e 



-t 



cos t sin t 
— sin t cos t 



P3.ll The state transition matrix is 

(2 1 - l)e -t —2te~ t 

m = 

2te~ t (-2i + l)e _t 
So, when xi(0) = x 2 (0) = 10, we have 

x(t) = <b(f)x(0) 
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or 



x\ (t) = lOe t 
X 2 (t) = 10e _t 

P3.12 (a) The phase variable representation is 



0 


1 


0 




0 


0 


0 


1 


x + 


0 


-48 


-44 


-12 




1 



y = [ 40 8 0]x . 

(b) The canonical representation is 



-2 


0 


0 




-0.5728 


0 


-4 


0 


z + 


4.1307 


0 


0 


-6 




4.5638 



y = [-5.2372 - 0.4842 - 0.2191]z 
(c) The state transition matrix is 



m 






where 





e" 6t - 3e _4t + 3e _2t 




|e- w - 2e~ 4t + |e _2t 


$iC 0 = 


—6e~ et + 12e _4t - 6e~ 2t 


^2(t) = 


-\e~ u + 8e- 4 * - |e- 2t 




36e~ et - 48e _4t + 12e“ 2t 




27e~ 6t - 32e _4t + he~ 2t 



$3 (t) 



Ip-6 1 _ 1„— 4 1 
8 e 4 e 

_| e -6t + e -4t 
| e -6t _ 4e -4 t 



+ h~ 2t 

_ Ip— 2t 
4 e 
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